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Abstract 

We introduce a momentum dependent mixing angle a(p^) which allow us to 
diagonalize at any external momentum p the one-loop corrected inverse propa- 
gator matrix of two coupled scalar fields while keeping the full momentum de- 
pendence in the self energies. We compare this method with more traditional 
techniques applied to the diagonalization of coupled scalars at the one-loop level. 
This method is applied to the CP-even Higgs sector of the Minimal Supersym- 
metric Model, defining the momentum dependent mixing angle a(p^), and cal- 
culating the two CP-even Higgs masses and the mixing angle at these two scales. 
We compare the results obtained in this way with alternative techniques. We 
make explicit the relation between a(p^) and the running mixing angle. We find 
differences between the mixing angle calculated with our method and the one 
calculated with more traditional methods, and these differences are relevant for 
Higgs searches at LEP2. 



1 Introduction 



In field theory it is common to find mixing between different species of scalars, fermions, 
or vector bosons. In the Standard Model (SM) we have the mixing between the gauge 
fields B and corresponding to the groups U{1) and SU{2) respectively. After a 
rotation given by the weak mixing angle ta.n6w = g'/g at tree level we find the mass 
eigenstates 7 and Z. Similarly, mixing between scalar particles are typical of two 
Higgs doublets models, and in supersymmetric theories mixing between fermions are 
also common (charginos, neutralinos) . 

In all cases it is trivial to find the mass eigenstates at tree level, however, one-loop 
radiative corrections will mix the tree-level diagonalized states (for example, Z and 7 
mix at one loop with charged particles in the loop). Moreover, the sum of the mixing 
graphs will be infinite and momentum dependent. In order to remove the infinities, 
mass counterterms and wave function renormalization constants are introduced. 

Several years ago, Capdequi-Peyranere and Talon studied the wave function 
renormalization of coupled systems of scalars, fermions, and vectors. Their approach 
include conventional mass counterterms and wave function renormalization plus a non- 
conventional field "rotation" that allow them to impose no mixing between the states 
at different scales. These scales are the masses of the different states. In this paper we 
generalize this idea and, at the same time, explain the nature of this field "rotation" 
by introducing a momentum dependent mixing angle between two coupled scalars. 

In the Minimal Supersymmetric Model (MSSM) the CP-even Higgs sector consists 
of two coupled scalars H and h, whose masses satisfy rriH > Tn^. Tree level mass 
relations are simple and specified by two unknown parameters: the mass of the CP-odd 
Higgs niA and the ratio of the two vacuum expectation value of the two Higgs doublets 
tan (3. Nevertheless, radiative corrections strongly modify the tree level mass relations, 
and since the lightest Higgs might be the first particle detected in the Higgs sector, it 
is necessary to achieve a good understanding of the effect of radiative corrections on 
this system of coupled scalars. For this reason, the CP-even sector of the MSSM is the 
best place to apply the method of diagonalization of coupled scalars we are proposing 
here. 

2 Renormalization of Coupled Scalars 

2.1 Conventional Wave Function Renormalization 

Similarly to ref. [|I| (the only difference is that they have m\2b = 0), consider the bare 
lagrangian corresponding to a system of two scalars: 

A = ^XibiP^ - mjf^)Xib + |x2fe(/ - mlf,)X2b - Xibmj2bX2b (1) 

If we denote by —iAfjlp^), i,j = 1,2, the sum of the one-loop Feynman graphs con- 
tributing to the two point functions and, after shifting the bare masses by m^^ —>■ 
mf — 6mf, i = 1,2, 12, and the fields by Xib (1 + \^Zi)xi, the effective lagrangian is 

^e// = i(Xl,X2)Sx(^^) (2) 
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where is the inverse propagator matrix, with matrix elements given by 

= -ml + {p^ -ml)6Zi + 6ml-A'[^{p^) 

^Up') = p'-ml + ip'-ml)SZ2 + Sml-AUp') (3) 

^Up') = -ml,-lml,{6Z, + 6Z2) + 6ml,-AUp') 

Although it is not a necessary assumption, in order to compare more easily with ref. 
we will assume for the moment that the two scalars are decoupled at tree level, i.e., 
= 0. In this case, if we want rrii, i = 1,2, to be the physical masses (the pole of 
the propagators) then the two mass counterterms are fixed through the relations 

Sml = AUml), 6ml = AUml) (4) 

Similarly, we may want to fix the wave function renormalization constants by setting 
to one the residue of each propagator at its pole. In this case we find 

6Z, = A'U^l), 6Z2 = A'Uml) (5) 

where the prime denote the derivative with respect to the argument. We may want to 
fix the Smi2 counterterm by imposing no mixing between xi and X2 at a given scale, 
for example at p^ = ml. In this case, the off diagonal element of the inverse propagator 
matrix is 

^Up') = Smj, - AUp') = AUml) - A^p') = -AUp') (6) 

At this point all the counterterms are fixed, and since Ai2(p^) is zero only at = ml, 
the two fields are not decoupled at a different scale. In particular, at the scale given 
by the mass of the second scalar, p"^ = m\, the mixing is not zero: 

T.Uml)=AUml)-AUml)^Q (7) 

Of course, this non-zero mixing is of one-loop order, and if the calculation of the one- 
loop scalar masses is done by diagonalizing the inverse propagator matrix order by 
order in perturbation theory, then this non-zero mixing would be a two-loop order 
effect. Nevertheless, this perturbative diagonalization can introduce large errors if 
radiative corrections are large. And this is the case with the CP-even Higgs sector of 
the Minimal Supersymmetric Model (MSSM). 

2.2 Mixed Wave Function Renormalization 

According to the previous section, the conventional wave function renormalization gives 
an inverse propagator matrix which is diagonal only at one particular scale. This is the 
motivation for the authors in ref. |1|] to define the following wave function renormaliza- 
tion xib (l-ai)Xi-AX2 and X2b ^ (1 -a2)X2 -/32Xi instead of Xii, {l + \5Zi)xi 
we use here. If we perform this transformation in the bare lagrangian in eq. (|l]) (and 
taking m^26 = in order to follow ref. [|l|) we find the following inverse propagator 
matrix elements 

E^i(/) = p^ - ml - 2{j)'^ - ml)ai + 5m\ - A\^(j)^) 

^Up') = p'-ml-2{p^-ml)a2 + 6ml-AUp'') (8) 

^Up') = -Pi{p'-ml)-(32{p'-ml)-AUp') 
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where the bar in Sis to differentiate with the matrix elements in eq. @. Setting to 
one the residue of the pole of each propagator they find 



1 ./X, 



1,2 



(9) 



where we have also included the relation between their and our 6Zi. Imposing no 
mixing between the two fields at p"^ = ml and also at p"^ = m\ they get 



^12 ("^ 



2 



2 ' 

mi 



^2 



AX (^2\ 



n2 



m. 



ml 



(10) 



This procedure is equivalent to take the inverse propagator matrix in eq. in the 
special case where m\2 = and 5m^2 — O5 perform a "rotation" (it is not a field 
rotation in the usual sense, it is just a wave function renormalization that mixes the 
two fields) in the following way 



1 P2 
Pi 1 



1 Pi 
P2 1 



:ii^ 



with the (3i being divergent, as it can be seen from eq. (|T0|). In this case, the inverse 
propagator matrix in eq. (H) is simultaneously diagonal at the two different scales 



p = mf and p 



Nevertheless, there will be a non-zero mixing between the two 



scalars at any other scale. 

It will be instructive to modify the calculation done in ref. [Q just presented in this 
section by considering — but 5m^2 7^ 0- this case, the only modification to the 
inverse propagator matrix is in the off-diagonal element in eq. (|]), and now is equal to 

^Up') = -Plip' - mj) - Piip' - ml) + 5ml, ' AUp') (12) 

The new coefficients /?/ can be calculated in the same way as before, that is imposing 
no mixing at the scale p^ = m\ and p^ = m^. We get 



P{ 



A^2 1 ^ 



2 2 
mf — m. 



Pi 



A\2{m\) - 5m\ 



12 



2 2 
m, — mi 



(13) 



The freedom introduced by the new counterterm 5m^2 gi^^ us the opportunity to cancel 
the divergency present in A^2{p'^)- This explains the superscript "/" in the constants 
p( : they are finite. In this way, the numerators in eq. (|13D are just the renormalized 
two point function Ai2(p^) evaluated at two different scales. The p( are then 



Pi 



mt, 



2 ' 
mi 



Pi 



AX (,^2\ 
■rl^ 



2 2 

m, — ml 



(14) 



This time, this wave function renormalization is equivalent to take the inverse propa- 
gator matrix in eq. (pi) and perform the following "rotation" : 



1 Pi 



[pi 



Pi 



Pi 
1 



(15) 



Again, the rotated inverse propagator matrix is diagonal only at the scales p^ 
and p^ = m\. 
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2.3 Momentum Dependent Mixing Angle 



In this paper we introduce a momentum dependent mixing angle a(p^) which will allow 
us to diagonalize the inverse propagator matrix at any momentum p|, ^ . Considering 
the already finite inverse propagator matrix elements in eq. (|]) (now we work on the 
general case rriu 7^ and 5m^2 7^ 0)) define the momentum dependent mixing angle 
by 



sin 2a (p^ 



cos 2a (p^ 



(16) 



ni{p')-n2{p')r+m2{p')Y 



The matrix TP^ijP') is diagonalized at any momentum by a rotation defined by the 
angle aijP') 



EX 



Ca(p^) 
Sa{p^) 



-SaiP^) 
Caip^) 



(17) 



-Sa(p^) Co,{p) 

where Sa and Cq, are sine and cosine of the momentum dependent mixing angle a(p^). 

In order to find the connection between the mixed wave function renormalization 
introduced in ref. |I| and the momentum dependent mixing angle introduced here, 
consider eqs. (^) and (16). In the case where = we find in first approximation 



Xp' 



aUp') 



ml 



m 



2 ' 



where (p^) is defined in eq. (P). Therefore, the mixed wave function renormalization 
defined by eq. ( P^ can be obtained from the rotation by an angle a(p^) defined in 
eq. (0) if we approximate eq. (|16]) to one-loop, order by order in perturbation theory, 
and evaluate the non-diagonal entries in the first rotation matrix in eq. (^) at two 
different scales: Sa at = ml and —Sa at = m^. 

The use of a momentum dependent mixing angle is an alternative to define a coun- 
terterm for this angle. In fact, the renormalization procedure is carried out in the 
unrotated basis and no mixing angle is defined at that level. Similarly, instead of 
renormalizing couplings of the rotated fields to other particles, we renormalize cou- 
plings of the unrotated fields to those particles and after that we rotate by an angle 
a(p^), where p"^ is the typical scale of the process, for example, p^ = mf if the rotated 
field Xi is on-shell. Usually, working out the radiative corrections in the unrotated 
basis implies one extra advantage, and that is the simplicity of the Feynman rules. In 
the following we will illustrate these ideas by renormalizing the CP-even neutral Higgs 
masses of the Minimal Supersymmetric Model (MSSM). 



3 The Minimal Supersymmetric Model 

The radiative corrections to the Higgs masses in the MSSM have been studied by many 
authors in the last few years, using three different techniques: the renormalization 
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group equation (RGE) method, the effective potential, and the diagramatic technique. 
It was estabhshed the convenience of the parametrization of the Higgs sector with the 
CP-odd Higgs mass m^, and the ratio of the vacuum expectation values of the two 
Higgs doublets tan/? = v^jvx. The radiative corrections to the charged Higgs mass 
were found to be small ^, P], growing as m^, unless there is an appreciable mixing 
in the squark mass matrix: in that case a term proportional to is non-negligible 
0. The corrections to the CP-even Higgs masses are large and grow as m^, and 
have profound consequences in the phenomenology of the Higgs sector @, H ^ [T0| . 
Two-loop corrections also have been calculated and shown to be important ||11||. 



The MSSM has two Higgs doublets p 



with Vi being the vacuum expectation value of the Higgs fields Hi. In the CP-even Higgs 
sector there are two coupled scalar fields (x?, X2)- The radiatively corrected propagators 
of the CP-even Higgs bosons depend on the two point functions Afj{p'^) = 1,2), 
where is the external momentum squared. In perturbative diagonalization of the 
one-loop mass matrix it is consistent to consider as a constant equal to the tree level 
mass matrix element. Nevertheless, it has been shown that for large values of the top 
quark mass, the perturbative diagonalization of the mass matrix is not reliable. We 
will study the effect of those approximations, namely the perturbative diagonalization 
of the mass matrix and the replacement of the external momentum in the self-energies 
by a constant, and compare them with the use of the momentum dependent mixing 
angle Q!(p^) that diagonalizes non-perturbatively the inverse propagator matrix of the 
CP-even Higgs sector, and we find numerically the pole of the propagators keeping the 
full momentum dependence. 

3.1 The model at tree level 

We start reviewing the neutral Higgs sector of the MSSM in the tree level approxima- 
tion. The CP-odd mass matrix given by: 

where is a soft symmetry breaking term, tp = tan (3 = V2/V1, and ti {i = 1,2) are 
the tree level tadpoles, whose expressions are 

t2 = mlHV2 - ml^vi + \{g^ + g'^)v2{vl - vl) . (21) 

Here, m^^ = mf + {i = 1,2), mf are soft supersymmetry breaking terms, and 
/i is the mass parameter in the superpotential. At tree level and at the minimum 
{ti = t2 = 0), this matrix is diagonalized with a rotation of an angle f3. The tree level 
CP-odd Higgs mass is: 

iOo = ^ . (22) 

S/3C/3 
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where the subscript "0" refers to the tree level. On the other hand, the mass matrix 
of the neutral CP-even Higgs bosons is 



Mi 



-\9WspCp-ml^ 
\gW4 + ^l2/ti3 + t2/v2 



(23) 



where g'z = g"^ + g'^ and 



jf + V2- From now on, we will use the following notation 



for its matrix elements: ( M^)jj = mf^. The tree level masses are obtained setting the 
tadpoles to zero and eliminating favor of m\ using eq. (P^). The answer is 



imjj^h)o = U^l + ml) ±yiml- m\fc% + {m\ + mlfs% , (24) 
with a tree level mixing angle 



(tan2a)o = f(tan2/?. 



(25) 

where we have omitted the subscript "0" from m\. Next, we calculate the radiative 
corrections to this mixing angle a, and introduce the momentum dependent mixing 
angle ctijP'). 



3.2 One-loop corrections to the CP-even Higgs masses 

To find the one loop corrections to the Higgs masses and mixing angle, we take the 
mass matrix for the CP-even fields and replace the bare quantities as indicated by: 



Hi 



X - 
m — 

Vi — 



:^ (1 + lSZ,)H,. 

■ A - 5A, 
m — 6m, 

Vi - 6vi, 



i = 1,2 
A = g,g\e... 

i = 1,2. 



(26) 



The one loop effective lagrangian in the CP-even Higgs sector (we have not rotated the 
original fields) has the form: 



(27) 



with the matrix elements given by 



S?i(/) = 
mp') = 
mp') = 
where we are allowed to set t,- 



P 



m 



(28) 



in mf^, given in eq. 



31). The mass counterterms 



are: 



6mii 
6mi2 



2r 2 I 2xf^l2\ I r,/ 2 2\ 2 2*^^/3 



2(mi 



6ti 
+ — 

mz)SfsCp— + — 



t 



V2 



(29) 



S/3Cf36ml - Sf3C(^6[ — —) - {m\ + ml)si3Cf3{ci - sl)—^ 
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and to fix them we adopt an on-shell scheme. Consider first the quadratic terms in the 
CP-odd Higgs sector (after rotating by an angle P): 



(30) 



where the relevant matrix element is 



S^-'Ip') = [p'-m\ + 6ml-AAA{p')\ZA (31) 
^ p"^ ~m\ + {p^ - m\)5ZA + 5m\ - Aaa{p^) ■ 

and the mass counterterm is given by 

W, = s(—)+^^- + <=l- (32) 

and the wave function renormalization constant is 

Za = s^Zi + cJZa SZa = s^SZh^ + cldZn^ . (33) 

To fix the mass counterterm 5m\ we adopt the following on-shell renormalization 
condition: 

ReE^^(m^) = ^ 5m\ = ReAAA(m^) (34) 

which means that the parameter m\ has been defined as the pole of the propagator. 

Similarly to the previous case, the mass counterterms for the Z and H^-boson are 
fixed in an on-shell scheme: 

Sml = ReAzz{rnl) , 6m^ = ReAww{m^) , (35) 

i.e., the parameters rriz and mw are the pole masses. 

Tadpoles are zero at tree level {ti = 0), and to fix their counterterms it is convenient 
to impose vanishing tadpoles at one-loop as well. If we call —iT\^^ to the sum of all 
Feynman diagrams (the 1-point irreducible Green's function) contributing to the one- 
loop tadpole, then the tadpole counterterm is equal to the sum of all one-loop tadpole 
graphs 

SU = T^l\p' = 0), 2 = 1,2. (36) 

Now we calculate the counterterm 6tp. We define the parameter tan (3 through the 
m-e^ and rriu^ masses |jl3[ [[ 

= . 2 2 2 \c2(3Awwimw) + Ae^eA^lj " ^^.^ei^T'lj] , (37) 

where the last two self energies receive no contributions from loops containing top and 
bottom quarks and squarks. 

In this way, with eqs. (0), (p^D, (pSj), (|36|), and (|37| ) we can calculate the mass 
counterterms in eq. (|29|) , which have been fixed in an on-shell scheme where m^, m^y, 
and rriA are pole masses, and the tadpoles are equal to zero at the one-loop level. 

-'^For alternative definitions of tan/? see for example ^ 
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The pole masses of the two neutral Higgs bosons are the zeros of the inverse propa- 
gator matrix YP^{p'^), whose matrix elements are defined in eq. (pS]), therefore the Higgs 
masses satisfy the following relation 



^22 



m,: 



yx ^™2^ 



'12 



m,: 



where we call 



P 



2 2 



2 2 
2 2 



+ 6mii 



h,H 



aUp') 

\x /„2^ 



(3^ 



%ip') 



(39) 



h, H are 



-m\sl- m\cl + Sm^2 - ^22 (P^ 
(m| + m^)s^C/3 + 5m^2 - ^12 (p^) 

From here we see that the wave function renormalization constants Zi, i 
canceled from the eigenvalue equation, implying the independence of the Higgs masses 
on the Zi. In ref. |jl6| the Higgs masses are calculated in this way, although without 
defining a momentum dependent mixing angle. 

The eigenvalues of the matrix E'^(p^) are the inverse propagators of the two Higgs 
bosons, and they are given by 



mp') = miip")- 

where i = H,h and the + 



^p')] ± yi^up') - ^up')? + m2ip')? (40) 



^22\J^ 7J 2 

sign correspond to the field h {H). These inverse 
propagators have a zero at the physical mass of the Higgs field: S/i(m^) = and 
Eii/(m|^) = and we solve these equations numerically. 

The wave function renormalization constants Zj, i = 1, 2, have not been calculated. 
We impose that the residue of the propagator of the CP-even Higgs bosons h and H 
are normalized to one. From eq. (|40| ) we get 



+ 



+ 



^22 



ml 



Zi 



^22 



m 



H) 



_d_ 
d_ 

Qp2 



mp': 



tUp') 



[m 



H) 



(41) 



Z2 Z\ 

from where the wave function renormalization constants can be calculated. 

The matrix E^(/r) is diagonalized at a particular momentum by a rotation 
defined by the angle oi{jP'\ This mixing angle satisfy 



tan 2a (p^ 



2t\2{j?)Z'l^zT 



t.Up')Zi-t.Up^)z, 



(42) 



In ref. it is used a wave function renormalization of the type in ref. [|T| to renormal- 
ize the CP-even Higgs sector of the MSSM, obtaining formulas analogous to eq. (p^), 
and then it is defined the angle a at two different scales. Those two angles are special 
cases of the momentum dependent mixing angle a(p^) in eq. (^2]). 

It is worth to mention that, since the residue of the pole of the propagator of the 
Higgs A is not normalized to unity, when renormalizing processes with external Higgs 
A it is necessary to multiply by the finite wave function normalization 

-1/2 



a/2 



drf^AKP 



-1/2 



g^2 



12 



Z1Z2 



^11^1 "T ^22^2 



(43) 
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which up to one-loop order reduces to Z 



1/2 



Finally, if we make the expansion Z, 



1/2 



1 + \Aj,^{m\ 



'l + 5Z,Y" 



as expected. 
1 + \5Zi. We get 



Sf2(/) 



P 



(44) 



And from here we have checked that all the divergences cancel in each matrix element 
in eq. (||). 



3.3 Perturbative limit 

At this point, it is instructive to make a perturbative expansion of tan2a(p^) defined 
in eq. (^). If we call At2a = tan 2a — t2ao5 ^"^^ keeping only one-loop terms, we get: 



At. 



2a 



(-200 



(m 



'2/3 



12 



+ i(^22 ~ ^11 )''^ 



2a 



(45) 



The variation of the angle a in eq. (|45|) is the finite one-loop correction to the tree 
level angle defined in eq. ([25|), calculated perturbatively. Now, to clarify even more the 
meaning of the angle a{p^) in eq. (|4^ ) and the perturbative one- loop correction to cto 
in eq. (|15|), consider the inverse propagator matrix S'^ with matrix elements given in 
eq. (|39D. If we rotate T?^ by an angle we will diagonalize the tree level part of the 
inverse propagator matrix, but not the one- loop part: 



R^mp')R, 



•tko 



P^ - - AhoHo{p^ 

-AhoHoiP^) 



-A 



ho Ho 



ip' 



P^ - ^lo - AhohoiP^) 



(46) 



where the rotation matrix is defined by 



COS ao sm 
- sin ao cos 



and m^^^ are the tree level masses given in eq. ( 
rotated CP-even fields: 



(47) 

ho and Hq are the tree level 

(48) 



and Aahip^) {a,b = ho^Ho) are the renormalized two-point functions. Now, to further 
diagonalize the inverse propagator matrix in eq. (^61) we need a rotation of one-loop 
order by an angle A a: 

1 Aft- 
-Aa 1 

Imposing that the off-diagonal matrix element in Rh.aRaQ'^^{p^)RZ^R^oi zero, we 
get 



(49) 



Aa 



C2ao^l2 + 1-3200(^22 — An 
m^Ho - ^ho 



(50) 
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where C2qo = cos2ao and S2q = sin2ao. Using the relations As a = Ca^Aa, and 
Acq = —SaoAa, we can prove that 

and from eq. (|50D and using S2ao("^li'o ~ ""^L) ~ '~i^A + ^'z)^w we recover eq. (^). 
In conclusion, we have proved that 

■Ra(p2) ~ RauRoo (52) 

2. e., a rotation by the angle ao followed by a rotation by the angle Aa is the perturbative 
approximation of a rotation by the angle C({p'^)- 



3.4 Relation with the running mixing angle 

Another useful point to clarify here is the relation between the momentum dependent 
mixing angle a(p^) and the running mixing angle a{Q), where Q is an arbitrary mass 
scale introduced by the momentum subtraction scheme DR in dimensional reduction 
[|17| 0. The renormalization group equation (RGE) satisfied by the angle a is directly 
related to the divergent terms of its counterterm. Since t2oo defined at tree level by 



t 



2m 



2ao 



X2 
12 



X2 X2 
"^11 - "^22 



its counterterm, defined by the relation t2ao = 

6t2a 2 



St2a, satisfy 
6m^2 + U^m^2 - 5mll)t2a 



(53) 



(54) 



^2q {'m\ + m%)s2i3 
This term is contained in the finite shift of the angle a given in eq. (|45|) , then we have 

At2a 5t2a 2 



2a 



t 



2a 



{m\ + m|)s2/3 



i(m^2'^^i/2 - m\^5ZH^)t2a + 



m\l{5ZH,+SZH,) + 



A 



hoHo 



(0) 



C2c 



(55) 



Since At2a is finite, to find the divergent terms of the counterterm St2a it is enough to 
look for the divergences of the second term in the right hand side of eq. (|55|) . In order 



to compare with ref . [T^ , where the contribution from the top quark to the counterterm 
for the angle a is calculated, we concentrate only in the top quark contribution. From 
eq. (^1]), and keeping only one-loop divergent terms, we have 



sz. 



Hi 



top 

div 



0, 



sz 



Hi 



top 
div 



A 



(56) 



where Nc 
given by 



3 is the number of colors, A is the regulator in dimensional regularization 

2 



A 



+ ln47r - 7^;, 



(57) 



^Or MS in the case of dimensional regularization used in non-supersymmetric theories. 
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n is the number of space-time dimensions, and 7^; is the Euler's constant. On the other 
hand, the contribution from the top quark to the mixing between and Hq is 



A 



(0) 



top _ 3Ncg^S2amf ^ 



div 



(58) 



Replacing eqs. (|^) and ( |58D into eq. (|55D, and considering that At2a is finite, we get 



St 



2a 



- top 



m 



- div 



(59) 



and this is the same answer we find in ref. |T^, with the exception of the sign, since 
here we define the counterterm of the angle a with the opposite sign. 

Now we have checked that we reproduce the correct divergent terms for the tan 2a 
counterterm, we turn to the relation itself between the renormalized tan2a(p^) and 
the running tan 2a (Q). In the MS scheme, each counterterm is fixed to cancel the 
divergent pieces of the corresponding loop corrections. In the case of tadpoles we have 



5tf ^ = [T,(g)]'^™ and - 5tf' + T,{Q) = T^'''{Q) . 
Therefore, the running tadpoles are equal toQ 



(60) 



ml^v, - ml,V2 + \{g' + g")vM - vl)\ (Q) = -ff ^(Q) , 
t^iQ) ^ [ml^V2-ml,v^ + l{g' + g'')v2{vl-vf)\iQ) = -fi^{Q). (61) 

Considering the tree level neutral CP-even Higgs mass matrix in eq. (p3|), we find that 



the renormalized inverse propagator matrix in the MS scheme has the following matrix 
elements: 



1 



msi{Q)+-friQ)-Air{p^Q) 

P' - mz{Q)sl{Q) - m^y^j^i^y 
[m|(Q) + ml{Q)]sp{Q)cp{Q) - I^r(/, Q) 



.2 ^2.^^„2.^^ -\(Q)4(g) + -TF(Q) - I^F(^)^ Q) (62) 

V2 



and from here we can find the renormalized a{p ) in terms of MS quatities 



tan2a(p^) 



(63) 



K(Q) + mims2m - '^AI^''{P\ Q) 



K(g) - mimc2p{Q) - I^^(p^ q) + ir^(p^ q) + irfr{Q) - i-fr'{Q) 

or, in first approximation 

2 



tan2a;(j9 



t2a{Q) 



{m\ - m|)c2/3 



At''{p\Q) 



(64) 



'ai^'{p\ q) - A\r{p\ Q) - -f^^^Q) + -rr (Q) 

V2 Vi 



'-2a 



To see the effect the running tadpoles have in the determination of running parameters, see mM 
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where the running tan 2a is defined as 



The relation between tan2a(p^) and tan2Q;(Q) in eq. 
the on-shell definition of the Z and A masses: 



(65) 

is completed when we give 



m'z = + ReAflKz, Q), = ml(Q) + Re<i(ml, Q), (66) 

where the tilde over the self energies means that the counterterms have already canceled 
the divergent terms. The relation between the on-shell definition of t/3 and the running 
tf^^Q) can be deduced from ref. . 



3.5 Numerical Results 

In this section we compare numerically the momentum dependent mixing angle method 
of diagonalizing coupled scalars with more conventional methods. We have chosen to 
compare with: 

(a) the tree level approximation; 

(b) the approximation where we diagonalize the matrix formed by the tree level quan- 
tities plus the pieces of the radiative corrections proportional to |2y, |^. By doing 
this we obtain the eigenvalues 



i(mi + m| + At) ± yiiml - m\)c2p - At? + {m\ + mlfslp , (67) 
and mixing angle 

tan2a= , ^fA + rnl)s,, 

{m\ - m|)c2/3 + At 



where: 



Sg'^mj 



1 1 
mi mf 

In ; (69) 



and 

(c) the leading logarithms approximation. 

The exact formulae for loops involving top and bottom quarks and squarks are given 
in the appendix (see also ref. p[), and loops corresponding to the gauge bosons, Higgs 
bosons, charginos and neutralinos are treated at the leading logarithm approximation. 
It is useful to evaluate the formulae in the appendix in the limit Mg^gy ^ ^? ^ "^i ~ 
^ ml and fi = Au = A^ = 0, where all the squark soft supersymmetry breaking 
mass parameters are assumed to be of the order of MgusY- The leading logarithms 



obtained in this way are in agreement with ref. [21 



— ^(i^t + n) In —2 + U{g mz) 

OZTT Cy^ '"'weak 
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vX / 2\ 2 2 2 2 2 3(7^ f'^'^t 2 2 1 i ^SUSY 

iOTT m^y Y / ^weak 

'J^'' "-Vy '"'weak 

where Pt = 1 — 4etS^ + 8e^s^ and = 1 + 4e6S^ + 8e^s^. Note that it also displayed 
the largest of the non-leading logarithm terms; it is proportional to the momentum p"^ 
and to the second power of the top quark mass. 

The hZZ coupling in the MSSM relative to the same coupling in the SM is given 
by the parameter sin(/? — a). In Fig. 1 we plot this parameter as a function of tan/?. 
The upper dotdashed curve correspond to the tree level approximation. The A^- 
improved approximation is in the dotted line [calculated with a defined in eq. 
The leading logarithms approximation is plotted in the lower dotdashed line. The 
parameter sin(/5 — a) calculated with the momentum dependent mixing angle a(p^) is 
plotted in the dashed curves for = m\j and in the solid curves for = ml, for the 
cases: {&) A = fi = 1 TeV and (b) A = —fi = 1 TeV, which define the squark mixing. 

From Fig. 1 we can learn that the tree level approximation can give completely 
wrong results. The A^-improved and the leading logarithms approximations give re- 
sults that are quite close to each other, indicating that the mf terms are the main 
contributions to the leading logarithms for the chosen parameters of this figure. Both 
approximations can differ from the results using a(p^), specially in the high tan/3 re- 
gion. The parameter calculated with our method shows a strong dependence on the 
squark mixing parameters, as it can be appreciated from the solid and dashed curves in 
cases (a) and (b). We also can see that the angle Q;(p^) has a small variation between 
the two physical external momenta = ml and p"^ = m\j. 

The same kind of graph is in Fig. 2, where we plot sin(/3 — a) as a function of 
tan (3. As opposed to the previous figure, here we consider lighter squarks: Mq = 
Mu = Md = MsusY = 200 GeV, Au = Ad = A = 140 GeV, and /i = -70 GeV. 
The most interesting feature here is that the parameter sin(/? — a) is substantially 
different in the two different external momenta p"^ = m\ and p^ = m'jj when tan (3 is 
large. In addition, this time not only the tree level calculation gives wrong results, 
but also the leading logarithms approximation. The fact that squarks are light implies 
that leading logarithms of particles other than squarks are also important, and this 
is reflected in the fact that the A^-improved approximation is very different from the 
complete leading logarithm approximation. For the parameters chosen in this figure, 
we see that the A^-improved curve is close to the curves calculated with a(p^), but 
this is an accident as we can see in the next figure. 

In Fig. 3 we plot sin(/5— a) as a function of for a fixed value of tan /? = 40 and for 
light squarks. Here it become evident that any of the three traditional approximations 
can give a value of sin(/? — a) with an error of 40% or more. Considering that the 
relevant parameter for the Higgs search is the MSSM is sin^(/? — a) we see that the 
error on the cross section can be of the order of 60% ! 
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The Higgs mass calculated with our method was compared to more traditional 
methods in ref . . In that reference we mention the differences between the Effective 
Potential method, the Renormalization Group Equations method, and the Diagramatic 
method. Here we discuss quantitatively these differences through an example given by 
the choice of parameters in Fig. 4. 

The Higgs masses are calculated by finding the zeros of the determinant of the 
inverse propagator matrix H^dr ) as indicated by eq. (0). The non-trivial momentum 

dependence of the matrix elements is in the self energies Usually, 

the momentum squared in the self energies is replaced by a constant. For example, in 
the effective potential method the momentum squared in the self energies is replaced by 
zero. Another typical choice is to take equal to the tree level mass. In Fig. 4 we have 
replaced by the same constant in the three self energies and plot as a function 
of this constant the two Higgs masses rrih and rriH calculated after that replacement. 
In the solid (dashed) line we have rrih {rriH) as a function of the squared root of the 
argument of the self energies. 

We see that, for the choice of parameters in Fig. 4, the dependence of the Higgs 
masses and tjih on is quite strong, specially for large values of p^. The dotted line 
is the diagonal defined by m = ^J^- The intersection of this line with the solid and the 
dashed curves give us the pole masses rrih and mn calculated with our method. These 
values are rrih = 64.8 GeV and rriH = 95.4 GeV. On the other hand, the intersection 
of the solid and dashed curves with the vertical line defined by p = corresponds 
to the masses calculated with the second derivative of the effective potential. These 
values are m'f/^ = 60.8 GeV and m^-^ = 108.4 GeV, and they are different from the 
pole masses. Therefore, it is clear that if we replace the momentum in the self 
energies by a constant, the Higgs masses calculated in that way may depend strongly 
on that choice. To be complete, we give the value of the Higgs masses calculated in 
(a) the tree level approximation m^^*^*^ = 90.7 GeV and m^*^*^ = 100.4 GeV, (b) the 
Aj-improved approximation m^* = 99.9 GeV and m^* = 129.3 GeV, and (c) the 
leading logarithms approximation m^f/" = 100.0 GeV and m^' = 133.9 GeV, valid for 
the choice of parameters in Fig. 4. The reason why the Higgs masses calculated with 
these approximations differs so much from the masses calculated with our method is the 
large value of the squark mixing: approximations (a) to (c) do not treat appropriately 
the squark mixing. An approximated formula which treats the squark mixing can be 
found for example in ref. p2[. As it can be seen above, a better value is obtained 



with the effective potential, but still, differences can be of 6% to 14%. It could be 
argued that the effective potential method can be improved by correcting the fact 
that it is found at zero external momentum, i.e., corrections of the type Am^ = 
A.hhijn'h^tree) ~ ^hhiS^) ■ Nevertheless, this correction implies that all the information 
comming from the effective potential is canceled out, and we are left with the pure 
diagramatic method. Therefore, we could have started with the diagramatic method 
in the first place and forget about the effective potential. The necessity of the above 
correction disappears in first approximation if m^'"'^^ = 0, as it was done in ref. (see 
also discussion in ref. II2F 
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4 Conclusions 



We have developed a new method of diagonahzing two coupled scalars by introducing 
a momentum dependent mixing angle a(p^), where p is the external momentum of 
the two-point functions. The dependence of the mixing angle on the momentum 
indicate us that the rotation matrix which diagonalizes the inverse propagator matrix is 
different whether we are at the pole mass = ml or = m^j or at any other scale. We 
have compared this method with the conventional wave function renormalization and 
with the mixed wave function renormalization In fact, we introduced the momentum 
dependent mixing angle as a generalization of the previous methods. 

We applied this method to the diagonalization of the CP-even Higgs bosons inverse 
propagator matrix in the MSSM. We used the diagramatic method in an on-shell 
renormalization scheme, where the tadpoles are exactly zero at one-loop, the masses 
m^, mvi/, and are defined as the pole masses, and tan/5 is defined through the 
on-shell definition of the slepton masses m^^ and m^;^. We calculate the wave function 
renormalization constants Zi and Z2 by imposing that the residue of the propagators 
of the CP-even Higgs bosons are exactly one. We do this by using a formula for Zi 
and Z2 valid for any number of loops and, therefore, specially useful when radiative 
corrections are large. We make explicit the relation between the momentum dependent 
mixing angle, which includes some effects of higher order loops, to the mixing angle 
calculated in the exact one-loop perturbative limit. We also make explicit the relation 
between the momentum dependent mixing angle a{p^) and the running mixing angle 
a;(Q), where Q is the arbitrary mass scale of the MS scheme. We give some numerical 
results by calculating the parameter sin(/3 — a), which is the ratio between the ZZh 
coupling in the MSSM to the same coupling in the SM. We compare this parameter 
calculated with the momentum dependent mixing angle a(p^) with (a) the tree level 
approximation, (b) the Aj-improved approximation (including only terms proportional 
to mf), and (c) the leading logarithms approximation. We find important numerical 
differences between the different methods, and they are relevant for the Higgs searches 
at LEP2 in the region of parameter space where ttia = 0{mz)- 

Finally, we calculate the pole masses of the CP-even Higgs bosons rrih and ttlh- We 
show that they are exactly independent of the wave function renormalization constants 
Zi and Z2 only if we use the formulas in eq. (pl). The Higgs masses are calculated by 
finding the zeros of the determinant of the inverse propagator matrix [eq. (^)], and this 
is the direct consequence of the definition of the momentum dependent mixing angle 
We compare with the Higgs masses calculated in the three approximations 
described in the above paragraph and we find that for some choices of parameter space 
there are non-negligible differences between them, and therefore, in those cases our 
method should be used. 

It is worth to mention also that our method can be easily generalized to the diagonal- 
ization of more than two coupled scalars. Also the generalization to the diagonalization 
of coupled fermions or coupled vector bosons is straight forward. 
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Appendix 



In this appendix we display the exact one-loop formulae we need to compute the inverse 
propagator matrix for the CP-even Higgs fields, for loops involving top and bottom 
quarks and squarks. We exhibit separately the contribution from the quarks and from 
the squarks loops. 

To obtain the mass counterterms Smfj given in eq. (|29|), the mass counterterm of 
the Z gauge boson is required. It is given by 5m| = Azzi.^'z) where the contribution 
from top and bottom quarks is: 



tb 



327r2c2, 



Ztt\ 



Zbb\ 



(71) 



where we use the notation -B^** = B22{'mz',rnf,m'f), and similarly for the other Velt- 



t2 . 2 (i.-y^S'^ 

man's functions. The contribution to the Z self energy from top and bottom squarks 
is: 



Azziml] 



tb 



2^2 



W 
1 ^2 



+ 
+ 



-\sl + e,sin2Bit'^ - 4^) + {\sl + e,slf{2B^t'' - A'^) 
■slcl{ABlt'' - 4 - 4) + \slcl{ABlt'' - 4 - K')] ■ (72) 



Here we use the notation st = sina^ and q = cos at where at is the rotation angle 
necessary to diagonalized the top squark mass matrix. Similar expressions are used for 
the sbottom mixing angle ab. 



The M^-boson self energy is 



Aww{^w) 



tb 



327r2 

for the quarks contribution, and 

tb NcQ"^ r 



AB. 



Wtb 
22 



- Al - 4 + (m^ - m 



2\ nWtb 
^b)^0 



Awwi''TT'w) 



327r2 



cictiAB^'^^^ - Al^ - A^^) + ci4i^B^'''^ - Al^ 



:>Wtib2 



1*1 



+ s'.cliABY^'^'^ - Al^ - aI^) + slsliABZ'-^- - 4 - A^^) 



(73) 

A'o') 
(74) 
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for the squark contributions, xxx The third mass counterterm we need to compute 
is 5{m\2/ SjjCp) and it is related to the CP-odd self energy and the tadpoles through 
eqs. (|3^) and (pif). Thus we need the following quantity: 



AAA{m\) 



.6ti 



6U 



-itb 



for top and bottom quark loops, and: 



/"^i -nAit I 2/1 T^Abb 



(75) 



AAAiniA) - sp— - Cp— 



327r2m2- 



w 



+ ml{fi + ADtpfBf^''^ + m,{^i/tp - ADtl)s^c,{Al^ - A^^)] (76) 

for top and bottom squarks loops. 

Finally, we need the two point functions A^y In the case of the self energies, we 
add the appropriate tadpole as indicated bellow. The contribution due to the top and 
bottom quarks is very simple: 



aUp') - - 
su 



aUp' 



V2 



aUp'' 



tb 



tb 



tb 



Ncg'^mf 
0. 



(77) 



where we notice that Xi couples only to top quarks, and X2 couples only to bottom 
quarks, implying that the mixing is zero. 

The contribution due to top and bottom squarks to the XiXi self energy minus the 
Xi tadpole is: 



aUp'' 



6ti 



tb 



2sp 



S2t 



+ 



+ 



+ 



+ 2 



1 2 



^zCpKo^t + (itSwC2t) + S2t 

2/1 2 \ I '^tfJ' 
"^zC/3(4 - etS^,)S2t + C2t 



^pt2t2 



Bt'^ + ^^2.(4^ - 4^) 



^\cp{\cl + ebs'^C2b) - 



^zCpi^sl - ebS^C2b) + 



mbAo ml 

— S2b 

2C/3 0/3 



rUbAD 



n 2 



Bo 



+ 2 



m 



Icpij + ebs'^)s2b + """^^ 



2cf, 
nibAi 



-S2b 



mt 



C/3 



B, 



pb'ib'z 



2c, 



C2b 



5, 



pb\b2 



rribA 



(7J 



D 



4.2 



S2b{Al^-Al^) 



where we use the notation S2t = sin(2aj), C2t = cos(2at), and similarly for a^. Similarly, 
the contribution to the X2X2 self energy minus the X2 tadpole, due to top and bottom 
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squarks is: 



V2 



tb 



167r2 



m 



w 



1^ S2t H 

2s/3 



5; 



+ 



2 / 1 2 , 2 N , "^t^t/ 
^zSl3\2h + etSwC2t) + — S2t 



+ 2 



m 



?S/3(-i + etslr)s2t 



mtA 



u 



2s 







-C2t 



ptit2 _ "^t^U 

4s| 



S2t{4-4) 



+ 



/9 



TDpbibi 
^0 



(79) 



+ 2 



'^\si3{-\sl + ebs'^C2b) - ^^S26 

2 / 1 , 2 \ , ^'-fo/^ 
^Z^m + ebSw)S2b + t; C26 



j^pb2b2 



1 2 



5, 







Finally, the X1X2 mixing due to squarks is given by: 



167r2 



m 



w 

X 



^zC/3(-2Ct + etSwC2t) + -7^ S2t 

,2 „ / 1 2 I 2 „ \ , ^tAu ^ , ""t 



fT^zSp^-^^t + etSvyC2t) + 



S2t + 

2S/3 5/3 J 



TDPtltl 

-Do 



+ 2 



2/1 2 \ I ^tA'" 

zC/^li - etS^)S2t + ^- — C2t 

/ 1 2 2 \ '^^^tf^ 

Z^l}\2^t + etSpyC2t) + -7^821 



^\sfi{-\ + etS^)s2t - C2t 

ZSfi 



ryptit2 
-^0 



2s. 



X 



rritAu ml 

— S2t 

2S/3 



Bq 



mlcpdcl + ebSwC2b) - 



rribAD ml 

—Z S2b 

2cp cp 



X 



"^Z^ 13(14 + ebSwC2b) - ^^S2b 



2c 



- 2 



2 /I I 2 \ , n^'b^D 
"n^Z^PKA + ^bSw)S2b + — C2b 



2c 



TDpbibi 
-^0 



2 /l I 2 \ I '^bA'' 
"^zS/^il + ebSw)S2b + C2b 



r>pb\b2 
^0 



+ 



2/12 2 N , ^bAp m^ 

n^zCpio^b - ^bSwC2b) + —p, S2b 

ZCi3 Cp 



X 



T^\si3{-\sl + ebs'^C2b) - ^^S26 



j^pb2b2 



(80) 



And this completes the formulas we need to calculate the momentum dependent mixing 
angle and the Higgs masses mh and mn in the approximation where only top 

and bottom quarks and squarks are considered in the loops. 
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Figure Captions 

Figure 1: Coupling hZZ relative to the SM coupling HsmZZ, sin(/5 — a), as a 
function of tan/5. It is plotted the tree level (upper dot-dashes), the At-improved tree 
level (dots), the leading logarithms (lower dot-dashes), and the parameter calculated 
with the momentum dependent mixing angle a(p^). In the last case, we use the mixing 
angle evaluated at the heavy Higgs mass scale a(m/^) (dashes) and the light Higgs 
mass scale a{mh) (solid), for two different choices of squark mixing: (a) A = fi = 1 
TeV and (b) A = -/i = 1 TeV. 

Figure 2: Couphng hZZ relative to the SM coupling HsmZZ, sin(/5 — a), as a 
function of tan/5. It is plotted the tree level (upper dot-dashes), the A^-improved tree 
level (dots), the leading logarithms (lower dot-dashes), and the parameter calculated 
with the momentum dependent mixing angle a(p^). In the last case, we use the mixing 
angle evaluated at the heavy Higgs mass scale a{mH) (dashes) and the light Higgs 
mass scale a{mh) (solid). We take tua = 100 GeV, Msusy = 200 GeV, A = 140 GeV, 
and /i = — 70 GeV. 

Figure 3: Coupling hZZ relative to the SM coupling HsmZZ, sin(/5 — a), as a 
function of m^. It is plotted the tree level (upper dot-dashes), the Aj-improved tree 
level (dots), the leading logarithms (lower dot-dashes), and the parameter calculated 
with the momentum dependent mixing angle a(p^). In the last case, we use the mixing 
angle evaluated at the heavy Higgs mass scale a(m//) (dashes) and the light Higgs 
mass scale a(m/j) (solid), for a squark mixing given by A = fi = 150 GeV. 

Figure 4: Neutral Higgs masses ruh (solid) and mn (dashes) as a function of the 
squared root of the argument of the self energies. We take tan/? = 50, tjia = 100 
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GeV, Mq = Mu = Md = Msusy = 1 TeV, Au = Ad = A = 1 TeV, and /i = -2 

TeV. The dotted line is the diagonal where m = \f^, and the intersection of this 
diagonal with the solid (dashed) curve gives us the pole mass m/j (mj^) in our approach. 
The intersection of the solid (dashed) curve with the vertical hne defined by = 
correspond to the Higgs mass nih {mH) calculated with the second derivative of the 
effective potential. 
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